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ABSTRACT:  The  electro-thermal  equations  describing  the  heating 
and  cooling  of  wire  bridge  electro-explosive  devices  are  solved 
for  constant- current  radar- type  input  pulses.  The  bridgewire 
temperature-time  history  is  obtained  for  a  variety  of  pulse 
amplitudes  and  repetition  frequencies.  Equilibrium  temperatures 
are  obtained  for  the  varied  input  conditions  and  are  all  com¬ 
bined  in  a  single  nomograph.  If  the  hot-spot  theory  is  assumed 
the  conditions  for  explosion  can  be  deduced.  For  instance,  for 
the  Mk  1  Scmlb,  a  temperature  rise  of  500°C  will  result  if  the 
value  of  I2rw  Is  1,000  microjoules  at  a  very  low  pulse  repeti¬ 
tion  frequency  (below  60  cps.)  or  is  100  microjoules  at  a  fre¬ 
quency  of  2,000  cps.  This  temperature  rise  would  be  expected 
to  cause  the  squib  to  fire  approximately  50$  of  the  time. 
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CHARACTERIZATION  OF  SQUIB, MK  1  MOD  0: 
Thermal  Stacking  from  Radar-Like  Pulses  (U) 


1.  INTRODUCTION 

1.1  The  mechanism  of  initiation  of  wire  bridge  electi’O- 
explosive  devices  (EEDs)  has  been  studied  at  this  Laboratory 
in  considerable  detail.  An  electro-thermal  model  has  been 
proposed  (reference  (l))  which  states  that  the  explosive  is 
triggered  when  the  wire  bridge  has  been  raised  to  a  particular 
temperature  by  "l2R",  or  ohmic,  heating  of  the  bridgewire . 
Experimental  data  obtained  to  date  have  not  required  an  altera¬ 
tion  of  the  basic  premise. 

1.2  These  studies  have  been  of  particular  interest  to, 

and  have  been  in  large  supported  by,  the  HERO  (Hazards  of  Electro¬ 
magnetic  Radiation  to  Ordnance)  program.  In  assessing  the 
vulnerability  of  an  EED  to  specific  electrical  environments,  it 
is  desirable  to  be  able  to  predict  the  probability  of  its 
response  to  various  types  of  electrical  inputs. 

1.3  The  present  report  deals  with  the  computation  of  the 
temperature-time  relationships  of  the  bridgewire  of  an  EED  when 
it  is  subjected  to  a  radar-like  train  of  RF  bursts  or  pulses. 

The  general  theory  of  the  electro-thermal  model,  the  digital 
computation  methods,  generalizations,  and  a  nomogram  of  typical 
results  will  be  given  in  this  report. 

1.4  The  Squib,  Mk  1  Mod  0,  has  been  studied  in  great  detail 
by  the  Naval  Ordnance  Laboratory.  It  is  considered  to  be  typi¬ 
cal  or  at  least  analogous  to  the  wire  bridge  EED3  in  general  use 
in  the  explosives  field.  Thus,  while  the  theory  and  detailed 
computations  have  been  applied  directly  to  this  particular  EED, 
it  is  firmly  believed  that  the  work  herein  can  be  applied 
directly,  or  with  little  change,  to  many  other  specific  EEDs. 

It  will  therefore  be  left  to  the  reader  to  make  such  adaptations 
or  extensions  as  may  be  necessary  to  suit  his  individual  needs 
rather  than  to  give  a  broader  and  more  comprehensive  treatment 
to  the  subject. 


1 


NOLTR  61-108 


2.  THEORY 

2.1  The  Squib,  Mk  1  Mod  0,  (shown  in  Figure  l)  has  been 
considered  as  a  lumped  parameter  system.  The  platlnum-iridlum 
bridgewlre  plus  its  immediate  surroundings  (explosive  when 
loaded,  air  before  loading,  or  explosive  simulant  when  dummy 
loaded)  has  a  heat  capacity  represented  by  Cp.  It  is  assumed 
that  Cp  is  independent  of  the  bridgewlre  temperature. 

2.2  The  dissipation  of  electrical  energy  in  the  bridgewlre 
causes  the  bridgewlre  to  be  elevated  in  temperature  by  an  amount, 
9.  The  bridgewlre  is  known  to  be  cooler  at  the  ends  than  in  the 
center  because  of  the  cooling  of  the  support  posts.  For 
simplicity,  the  bridgewlre  temperature  is  considered  to  be  an 
average  of  the  temperature  along  its  length.  The  hottest  point 
on  the  wire  will  be  the  most  likely  to  cause  Initiation.  Yet, 
experimentally,  this  hottest  point  does  not  differ  markedly  from 
the  average  temperature. 

2.3  A  steady  state  condition  of  electrical  energy  dissipa¬ 
tion  in  the  bridgewlre  will  lead  to  a  corresponding  steady  3tate 
flow  of  thermal  energy  into  the  support  posts  and  into  the 
explosive,  or  simulant,  in  contact  with  the  bridgewlre.  Under 
these  steady  state  conditions  the  heat  energy  flow,  and  therefore 
the  electrical  energy  input,  will  be  proportional  to  the  tempera¬ 
ture  elevation,  ©.  The  proportionality  constant,  y  ,  is 
designated  as  the  heat  loss  factor. 

2.4  These  parameters  can  be  combined  in  an  instantaneous 
energy-balance  equation 

Cp  _H_  +  r  9  =  P(t), 

where  P(t)  is  the  power  time  function. 

From  this  differential  equation  9  (t)  functions  can  be  derived 
which  describe  the  thermal  response  of  the  EED  bridgewlre  to 
various  electrical  input  wave-forms.  Experimental  verification 
of  the  equation  is  facilitated  by  the  technique  of  using  the 
EED  bridgewlre  as  its  own  resistance  thermometer.  The  change 
of  bridgewlre  resistance  with  temperature  can  be  expressed 
(using  again  the  idea  of  lumped  constants  and  of  9  as  the  average 
bridgewlre  temperature  elevation)  by 

R  -  R  (1  +  cL9) 
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where  R  Is  the  resistance  of  the  EED  at  a  stated  reference 
temperature . 

and  c L  13  the  corresponding  temperature  coefficient  of 

resistance  measured  for  a  temperature  change  from 
this  same  stated  reference  temperature. 

This  relationship  is  based  on  the  assumption  that  the  resistance 
varies  linearly  with  the  temperature  throughout  the  practical 
temperature  range . 

2.5  Because  of  the  low  resistance  of  the  EED  compared  to 
the  probable  magnitudes  of  source  impedances  that  might  be 
expected  in  weapon  circuitry  carrying  RF  power,  it  was  decided 
to  consider  the  case  for  constant-current  pulse  trains.  An 
additional  reason  for  this  choice  is  that,  for  EEDs  with  a 
positive  value  of  the  coefficient,  pC  ,  heating  of  the  bridge- 
wire  raises  its  resistance  and  therefore  increases  the  power 
it  absorbs,  even  though  the  pulse  amplitude  1b  unchanged.* 
Constant  current  pulse  data  would  therefore  provide  a  conser¬ 
vative  basis  for  making  system  safety  estimates.  The  differ¬ 
ential  equation  expressing  the  lumped  parameter  behavior  under 
constant  current  input  conditions  is: 

cn  +  re  -  i2r  (i  +*e)  (1) 

p  dt 


which  can  be  solved  to  yield 


where  W  is  the  pulse  length 

and  ©w  is  the  temperature  at  the  end  of  the  pulse. 

2.6  Whenever  the  bridgewire  temperature  is  above  ambient 
with  no  power  applied,  the  differential  equation  expressing  the 
bridgewire  cooling  will  be  found  to  apply: 

cD  _dQ  +  re  -  0.  (3) 

K  d  t 

♦Positive  values  of  dt  cause  constant  voltage  pulses  to  deliver 
less  power  as  the  temperature  rises,  i.e.. 
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Taking  time,  t,  measured  from  the  Instant  the  power  is  removed 
and  cooling  begins,  and  taking  the  bridgewire  temperature  to  be 
9m  at  the  instant  t  »  0,  equation  (3)  can  be  solved: 


6 


c 


9m  exp 


-  y  tc 

*  9m  exp 

4C 

cp  _ 

m  ~ 

r 

(4) 


where 

tc 

and  ©c 


n 

the  thermal  time  constant,  equals  p  , 

y 

is  the  cooling  time  between  current  pulses  , 
is  the  temperature  at  the  end  of  the  cooling  time. 


2.7  To  obtain  the  temperature -time  history  of  an  EED 
subjected  to  a  train  of  constant  current  pulses  it  is  necessary 
to  apply  equations  (2)  and  (4)  in  alternation  in  as  many  cycles 
as  may  be  appropriate  to  the  study. 


2.8  If  the  cooling  portion  of  the  cycle  1b  long  compared  to 
/,  the  bridgewire  will  return  to  ambient  temperature  at  the  end 
of  the  cycle.  In  this  case  each  cycle  will  repeat  the  history 
of  the  preceding  cycle.  If  the  pulse  repetition  frequency, 
that  is  the  frequency  at  which  the  pulses  are  repeated,  is 
increased  by  decreasing  the  cooling  portion  of  the  cycle  beyond 
a  certain  point,  then  at  the  beginning  of  the  pulse  train  the 
loss  of  heat  during  cooling  will  be  less  than  the  gain  during 
the  heating  portion  of  the  cycle.  In  this  case  the  temperature 
of  the  bridgewire  will  be  greater  for  one  cycle  than  it  was 
for  the  preceding  one.  As  a  result  of  this  Increase  in  temp¬ 
erature  both  the  rise  and  fall  of  temperature  during  the  heating 
and  cooling  portions  of  the  cycle  are  increased.  Howevp”,  sir-'e 
the  increase  in  the  rise  of  temperature  during  the  heating  portion 
of  one  cycle  with  respect  to  that  of  the  preceding  cycle  is 
linear  and  the  fall  in  temperature  during  the  cooling  portion 
Increases  exponentially,  the  loss  of  heat  due  to  cooling  will 
eventually  become  equal  to  the  gain  in  heat  during  the  heating 
portion  and  the  temperature  of  the  wire  again  goes  through  a 
stable  history  from  one  cycle  to  the  next  as  shown  in  Figures 
2  and  3.  If  the  input  energy  in  the  constant  current  pulses  is 
high  enough  then  the  bridgewire  temperature  will  rise  high  enough 
to  initiate  the  explosive  surrounding  it.  If  the  energy  is  not 
high  enough  the  bridge  temperature  will  stabilize  at  a  value 
lower  than  the  ignition  temperature  of  the  surrounding  explosive. 
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3.  COMPUTATIONS  AND  RESULTS 

3.1  An  IBM  704  Computer  was  programmed  to  express  the 
constant  current  heating  with  alternate  natural  cooling.  The 
Fortran  Program  is  given  in  Appendix  A.  The  input  data  for  the 
IBM  704  program  consisted  of  the  current,  the  initial  resistance, 
the  repetition  frequency,  the  pulse  width,  and  the  thermal  para¬ 
meters;  y  i  <k  >  and  Cp.  Using  these  data  and  equation  (2),  the 
program  computed  the  temperature  of  the  bridgewire  at  the  end 

of  the  heating  portion  of  the  cycle.  Having  obtained  this 
result  equation  (4)  was  used  to  find  the  temperature  at  the  end 
of  the  cooling  portion  of  the  cycle.  The  resistance  of  the 
bridgewire  which  corresponded  to  these  temperatures  was  also 
calculated.  The  program  also  computed  the  energy  absorbed  in 
the  heating  portion  of  each  cycle.  The  computation  was  then 
repeated  using  the  temperature  and  resistance  at  the  end  of  the 
cooling  portion  of  the  cycle  as  the  values  for  the  start  of  the 
next  cycle.  After  each  cycle  tests  were  made  of  the  maximum 
temperature  of  the  cycle,  the  increase  in  this  temperature  over 
that  of  the  preceding  cycle,  and  the  number  of  cycles  computed. 
These  values  were  compared  with  preassigned  values  of  these 
quantities.  Computation  ceased  after  the  preassigned  number  of 
cycles  or  if  the  temperature  exceeded  the  preassigned  value  for 
the  maximum  temperature,  which  usually  was  2,000°C,  or  if  the 
increase  in  temperature  over  the  preceding  cycle  was  less  than 
one  half  of  one  percent.  Otherwise,  computation  proceeded  to 
the  next  cycle. 

3.2  The  initial  resistance  in  almost  all  cases  was  taken  to 
be  1.0  ohm.  Pulse  repetition  frequencies  from  6u  to  3,000  cycles 
per  second  were  studied.  In  most  cases  the  on  time  for  the  pulse 
was  short  compared  with  the  total  time  of  the  cycle  and  small 
compared  to  the  thermal  time  constant  ^  .  Under  these  conditions 
the  cooling  portion  of  the  cycle  can  be  considered  to  be  independ¬ 
ent  of  the  pulse  width.  The  temperature  rise  during  a  pulse  will 
be  a  function  of  the  energy.  Thus,  for  any  given  pulse  frequency 
the  same  maximum  peak  temperature  will  be  reached  for  any  com¬ 
bination  of  current,  resistance,  and  pulse  width  for  which  the 
product  I^RW  Is  constant. 

3.3  Results  printed  out  (Figure  4)  included  the  resistance 
and  temperature  at  the  beginning  and  end  of  the  heating  portion 
of  the  cycle,  the  energy  of  the  heating  portion  of  the  cycle, 
and  the  total  energy  used  in  this  and  previous  cycles.  Compu¬ 
tations  were  carried  out  for  different  combinations  of  current, 
pulse  repetition  frequency,  and  pulse  width  selected  so  as  to 
give  maximum  peak  temperatures  of  2,000°C  or  less. 

3.4  It  was  found  that  if  the  logarithm  of  I  is  plotted  against 
the  logarithm  of  the  pulse  width  for  any  given  frequency  and 
(naximum  temperature,  the  result  is  very  nearly  a  straight  line. 

This  was  done  for  several  representative  fre q  u  e  ncie  s  and 
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is  shown  in  Figures  5  and  6.  In  each  case  lines  were  drawn 
representing  temperatures  of  300,  500,  and  1,000°C.  The  results 
of  the  calculations  can  all  be  combined  into  a  single  graph  in 
which  the  product  I2rw  is  the  ordinate  and  the  pulse  repetition 
frequency  is  the  abscissa.  The  results  can  be  shown  by  drawing 
contour  lines  representing  the  attained  maximum  peak  temperature 
(Figure  7).  In  order  to  locate  these  contour  lines  it  was 
convenient  to  use  a  method  of  computation  which  would  have  a 
given  temperature  as  part  of  the  input  data  and  calculate  the 
point  corresponding  to  this  temperature  for  some  given  combina¬ 
tion  of  the  other  variables. 

3.5  Appendix  B  gives  the  program  as  it  was  modified  to 
accomplish  this  result.  The  modified  program  was  arranged  to 
find  the  pulse  width  which  would  give  the  desired  maximum  peak 
temperature  for  a  given  pulse  repetition  frequency.  This  is 
done  by  repeated  interpolations  between  previously  determined 
pairs  of  values  of  the  pulse  width  and  temperature.  The  program 
begins  with  two  pulse  widths  and  their  corresponding  maximum 
temperatures.  It  then  interpolates  or  extrapolates  for  the 
pulse  width  corresponding  to  the  desired  temperature.  Using 
this  value  of  the  pulse  width  it  then  computes  the  maximum  peak 
temperature  as  in  the  previous  program.  This  temperature  is 
compared  with  the  desired  temperature.  If  the  difference  is  less 
than  5°C,  the  computation  ceases.  Otherwise,  the  new  pulse  width 
and  temperature  replace  one  of  the  previous  pairs  of  values  and 
the  process  is  repeated.  One  of  the  original  pairs  of  pulse 
width  and  temperature  values  can  be  taken  to  be  zeros;  the  other 
pair  may  be  merely  an  estimate  for  which  no  great  accuracy  is 
required.  The  process  will  find  a  satisfactory  value  for  the 
pulse  width  within  a  few  iterations. 

-  O 0 

3.6  Values  of  ICRW  for  frequencies  of  pulse  repetition 
frequency  from  twenty  to  two  thousand  cycles  per  second  giving 
maximum  peak  temperatures  of  300,  500,  and  1,000°C  above  ambient 
temperature  were  computed  in  this  way  and  are  given  in  Table  1. 
These  results  are  also  shown  graphically  with  contour  lines  for 
the  maximum  temperatures  drawn.  In  order  to  simplify  the  location 
of  the  ordinate  value  for  this  graph,  a  simple  nomogram  was 
designed  which  gjves  the  ordinate  for  known  values  of  I  and  W 
with  resistance  R.  The  nomogram  is  shown  in  Figure  8.  These 
computations  used  typical  values  of  the  thermal  parameters:  Y~  , 

600  microwatts  per  degree;  cL  ,  J05  parts  per  million  per  degree* 
and  Cp,  2.4  microjoules  per  degree.  A  few  calculations 

were  male  using  a  value  for  eC  of  900  parts  per  million  per 
degree.  Temperature  reached  under  these  conditions  are  higher 

*It  should  be  noted  that  d.  is  a  temperature  dependent  constant 
and  must  be  selected  to  correspond  to  the  base  temperature  from 
which  resistance  calculations  are  made.  If  the  base  temperature 
or  base  resistance  changes  then  oL  must  also  be  changed. 
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than  for  the  smaller  value  of  .  The  difference  is  greater 
when  comparing  the  maximum  peak  temperature  than  when  comparing 
the  first  pulses.  For  the  first  pulse  the  difference  is  of  the 
order  of  5^.  At  higher  frequencies  the  maximum  peak  temperatures 
differ  by  about  15$.  Utilization  of  the  computational  technique 
therefore  requires  a  realistic  choice  of  numerical  values  for 
the  various  parameters.  Experience  to  date  has  shown  that  the 
determination  of  n  is  less  subject  to  experimental  error  and  is 
less  variable  than  that  of  the  other  parameters.  The  determina¬ 
tion  of  y  and  Cp  is  the  least  precise  since  the  error  of  their 
determination  must  include  the  error  of  determination  of  . 

The  extremes  of  the  individual  values  of  the  most  variable 
parameter  rarely  exceed  0.6  to  1.5  times  the  central  value. 


4.  CONCLUSIONS 

4.1  The  results  obtained  in  this  investigation  have  been 
summarized  in  Table  1  and  as  a  nomograph  (Figure  8) .  The  table 
gives  initial  conditions  which  determine  the  pulse  size  required 
to  heat  one  particular  type  of  brldgewlre  ultimately  to  a  maximum 
peak  temperature  of  300,  500,  or  1,000°C  above  ambient  for  pulse 
repetition  frequencies  from  20  to  2,000  cycles  per  second.  The 
initial  conditions  are  most  conveniently  expressed  as  a  pseudo¬ 
energy  per  burst  term  in  units  of  microjoules  per  burst.  This 
term  is  computed  on  the  assumption  that  the  pulse  power  is  main¬ 
tained  at  its  initial  value  (the  power  level  before  any  bridge- 
wire  resistance  change  occurs) .  The  actual  pulse  energy  will, 

of  course,  be  greater  than  this  due  to  the  heating  and  corres¬ 
ponding  resistance  increase  of  the  brldgewlre. 

4.2  The  computed  results  can  be  used  to  assess  the  expected 
response  of  ar}  EED  to  a  particular  electromagnetic  energization. 

If  values  of  R  ,  oL ,  f ,  and  Cp  are  known  for  a  particular  EED, 
it  is  possible  to  compute  the  interrelationships  between;  (a) 
maximum  peak  temperature,  (b)  pulse  repetition  frequency,  and 
(c)  initial  pulse  magnitude  (which  in  turn  can  be  broken  down 

to  pulse  width  and  initial  pulse  amplitude).  For  instance, 
for  the  Squib,  Mk  1  Mod  0,  a  temperature  rise  of  500°C  would 
be  expected  at  a  pulse  repetition  frequency  of  2,000  cycles  per 
second  for  a  10-ampere  pulse  with  a  pulse  width  of  1.1  microsecond. 
The  same  temperature  would  be  attained  at  pulse  repetition  fre¬ 
quencies  below  approximately  30  cycles  per  second  for  the  same 
current  amplitude  provided  the  pulse  width  were  10.5  microseconds. 

4.3  Once  the  temperature  can  be  deduced  the  probability  of 
initiation  of  the  EED  can  in  turn  be  estimated  provided  informa¬ 
tion  is  available  as  to  the  temperature  which  will  cause  initia¬ 
tion  of  the  EED.  This  information  usually  can  be  found 
experimentally  from  DC  pulse  firing  data.  Inference  of  the 
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necessary  temperature  based  on  the  "Hot-Spot"  Theory  of  Bowden 
and  Yoffee  (reference  (a))  has  been  found  to  give  satisfactory 
results  (reference  (3)). 

4.4  Subsequent  to  the  work  herein  reported  a  technique  was 
devised  to  monitor  the  bridgewire  heating  and  cooling  of  EEDs 
being  energized  by  microwave  energy  from  a  9 -KMC  radar.  Quali¬ 
tative  and  preliminary  quantitative  comparison  show  agreement 
between  results  and  predicted  values  indicating  that  the 
phenomena  of  thermal  stacking  and  the  appearance  of  the  time- 
temperature  wave-forms  as  they  are  described  herein  are  in 
accord  with  experiment. 
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Table  1 

Relationship  Between  I2RW  (micro Joules ) , 

Pulse  Repetition  Frequency,  and  Equilibrium  Peak  Temperature 


PULSE 

REPETITION 

FREQUENCY 

(cps) 

EQUILIBRIUM  PEAK 

TEMPERATURE  ( 

°C) 

300 

500 

TODD 

20 

653.1 

1029.2 

1819.9 

40 

652.3 

1026.7 

1814.7 

60 

642.3 

1009.8 

1778.3 

100 

594.5 

928.4 

1608.4 

140 

530.9 

827.7 

1408.1 

200 

452.6 

695.1 

1156.0 

250 

398.3 

606.8 

1019.3 

300 

353.9 

539.4 

878.5 

600 

211.5 

315.7 

505.O 

1000 

137.7 

205.2 

324.3 

1400 

103.0 

152.6 

240.1 

2000 

74.6 

110.9 

172.4 
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NOTES: 

1.  IGNITION  BEAD  -  APPROX.  5  MG  DDNP/KCI03 

2.  FLASH  CHARGE  -  APPROX  45  MG  BLACK  POWDER 
3  BASE  CHARGE  -  APPROX  45  MG  BLACK  POWDER 
4.BRI0GE  WIRE- O.OOl"  PLATINUM-IRIDIUM  0.060"  LONG 


BRIDGE-WIRE  TEMPERATURE  (ABOVE  AMBIENT), 


10  20  30  40 

TIME,  MILLISECONDS 

NOTE:  CONSTANT  CURRENT  I4AMP, 

PULSE  REPETITION  FREQUENCY  499.5CPS 
PULSE  WIDTH  2  MICROSECONDS 

.2  THERMAL  STACKING  OF  THE  SQUIB,  MK  I  BRIDGE 


NOTE :  CONSTANT  CURRENT  PULSES  APPLIED  TO 
SQUIB,  MK  I  LOADED.  PULSE  REPETITION 
FREQUENCY  499.5  CPS  PULSE  WIDTH, 

2  MICROSECONDS. 


FIG.  3  THERMAL  STACKING  AS  A  FUNCTION  OF  PULSE  AMPLITUDE 
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APPENDIX  A 


A-l.  Fortran  II  Is  an  Automatic  Programming  System  of 
Formula  Translation.  It  is  a  method  of  rewriting  mathematical 
relations  in  a  symbolic  form  which  can  be  converted  by  a  high 
speed  digital  computer  into  a  productive  sequential  computa¬ 
tional  procedure.  Fortran  II  can  now  be  interpreted  by  a 
number  of  different  makes  of  high  speed  computers. 

A-2.  Pages  A— 3  and  A-4  are  a  photographically  reproduced 
copy  of  the  Fortran  II  listing  of  the  computational  program 
described  in  Paragraph  3.1.  An  exact  keypunch  duplication  of 
this  listing,  one  card  per  line,  indluding  all  spaces  within 
each  line,  should  produce  a  source  deck  which  can  be  compiled 
by  any  activity  to  whom  Fortran  II  is  available. 

A — 3 .  A  description  follows  of  each  of  the  program 
symbols  of  quantities  required  for  the  computations. 


Program 
Symbols : 


TMAX  Upper  bound  of  computed  temperature  of 

EED  brldgewire  in  degrees  Centigrade. 

PC  Ratio  of  maximum  temperature  attained 

in  one  pulse  to  that  of  preceding  pulse 
which  is  considered  to  represent  a 
stable  temperature  condition. 

NMAX  Maximum  number  of  pulses  to  be  computed. 

CUR  Constant  current  in  amperes. 

RZ  Initial  resistance  in  ohms. 

RATE  Pulse  repetition  frequency  in  cycles 

per  second. 

TP  Time  in  seconds  during  which  the 

constant  current  flows  in  each  pulse. 

GAM  Heat  loss  factor.  Gamma,  in  watts  per 

degree  Centigrade. 

ALF  Temperature  coefficient  of  resistivity. 

Alpha,  in  ohms  per  ohm  per  degree 
Centigrade . 

CP  Heat  capacity  in  Joules  per  degree 

Centigrade . 
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APPENDIX  A  (cont'd) 


Program 
Symbols : 


T1  Temperature  at  the  beginning  of  pulse. 

R1  Resistance  at  the  beginning  of  pulse. 

T2  Temperature  at  end  of  heating  portion 

of  cycle. 

R2  Resistance  at  end  of  heating  portion 

of  cycle. 

T3  Temperature  at  end  of  cooling  portion 

of  cycle. 

R3  Resistance  at  end  of  cooling  portion 

of  cycle. 

TLST  Maximum  temperature  of  preceding  pulse. 

TS  Length  of  cooling  portion  of  cycle. 

PENG  Energy  in  the  pulse. 

ENG  Total  Energy. 

TIME  Total  time  elapsed. 
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C  PROGRAM  FCR  FINCING  T EHP ER ATURE- T I  ME  HISTORY  OF  AN  EEC  SUBJECTED 

C  TO  CONSTANT  CURRENT  RACAR-LIKE  PULSES 

5  READ  INPUT  TAPE  0,301,  TMAX.PC.NMAX 
7  NN  =  C 

10  RE AC  INPUT  TAPE  C , 30 2 , CUR , RZ , R AT E , T P , GAM , AL F , C P 

11  IF (RATE  )  115,115,12 

12  IF t  l.0-TP*RATE  )  115,115,16 

16  T  1  =  C  .  0 

17  R  1  =  R  Z 

19  TL  S T  =  T  1 

20  N  =  0 

21  ENG'O.C 

23  TS=  1  . O/RATE- TP 

24  EXP  =  0.36787944»*( TS«GAM/CP  ) 

28  WRITE  OUTPUT  TAPE  0,304,  CUR , R Z , R A T E  ,  TP , G AM , AL F , C P 

29  NN  =NN+  5 

30  CS0=CUR*«2 

31  ENZ=RZ»CSO 

40  EN  =  R  1  »CSC 

41  D=EN»ALF-GAM 

42  C  =  C  P  /  0 

43  E*2.71828183«« (TP/Cl 

44  A=EN/D 

45  N  =  N ♦  1 

46  B  =  A • AL  F 
48  EC  =  E- 1 .0 

50  T2  =  T  1+EC»A 

51  T3=T2«EXP 

52  R2  =  RZ  * (  1 .0  +  ALF-T2  ) 

53  Ri  =  RZ» (  1.0  +  ALF»T31 

60  PENG=EN«TP+B«(C*EC-TP)*ENZ 
62  ENG=ENG*PENG 

74  AN=N 

75  T  I M  E  =  AN / R A  T  E 

76  I  F ( NN-40  )  79,79,77 

77  WRITE  OUTPUT  TAPE  0,307,  N , PENG 
NN  =  C 

78  GO  TO  80 

79  WRITE  OUTPUT  TAPE  0,308,  N , PENG 

80  WRITE  OUTPUT  TAPE  C  ,  309 , T I  ME , ENG 

81  WRITE  OUTPUT  TAPE  0,310 

e2  WRITE  OUTPUT  TAPE  0,311,T1,R1 

83  WRITE  OUTPUT  TAPE  0,3!2,T2,R2 

84  PX=(T2+273.CI/(TLST*273.0) 

86  IF1TMAX-T2)  7,7,87 

87  IF(PC-PX)  88,7,7 

88  IF(NMAX-N)  7,7,90 

90  R 1 =  R 3 

91  TLST=T2 

92  T 1  =  T  3 

93  NN=NN+6 
95  GC  TO  40 

115  CALL  RETURN 

1 1 6  STOP 

301  FORMAT ( F9. C , F9 .4 ,  I  9 > 

302  FORMAT ( 3F8.C,4E 1 1.4) 

304  FORMAT (  18  1  64  (  18  )  188T8ERMAL  P A R AM E T E R S/ 1  0 (  1  H  )  7  8HCURR  ENT  RESIST 

1ANCE  PULSES  PULSE  WICTH  GAMMA  ALPHA  CEE  PEE/1211  304 

2H  178HAMPS  C8MS  PER  SFC  MICROSECONDS  MICWTS/OEG  P/  304 
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3M/UEG  f'ICJLS/CEG/F  16.2.F  13.3  ,F  1  1  .2,6PF  14.3.2F  11  .  1  ,F  12.3) 

307  FCRP  AT (  1  h  1  1 6 (  1 H  )  1  lFCYCLE  COUN T I  1 C , 1 3 ( 1 H  HSHENERCY  IN  PULSE3PFI1 
13, 12H  PILLI JCLLES  ) 

308  FC  R  N AT (  1  HO  1  6 (  1 H  IIIFCYCLE  CCUN T  I  1  C ,  1  3 (1 H  )  1  bHENERG Y  IN  PULSE3PF11 
13, 12H  PILLI  JCLLES  ) 

3C9  FCRRAT  (  1  7  (  1H  )  12FELAPSEC  TIRE3PF  1  2 . 2 , 2bH  MILLISEC  TOTAL  ENERGY 
1F1  1 .3,  12H  f'lLLI  JCULES) 

3 1 C  FORMAT (  1 7 (  Ih  )  1  1 FT EH  PER A T URE2 3 (  1 H  ) 1 OHRE S  I  ST ANCE  ) 

311  FCRPAT(19(1F  )  12HBEFCRE  PUL SE F  1 C . 2 , 24H  DEG  C  BEFORE  PULSE 

1F12.3.5H  OHf'S) 

312  FCRKAT (  191  1H  111FAFTER  PU LSE F 11  . 2 , 2 3H  DEG  C  AFTER  PULSE 

1F13.3.5H  OHPS) 
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APPENDIX  B 


B-l .  Appendix  B  contains  the  Fortran  II  listing  for  the 
computational  program  described  in  Paragraph  3.5  together  with 
a  listing  of  symbols  used  in  this  program  which  were  not  used 
in  the  first  program.  The  list  of  symbols  given  in  Appendix  A 
applies  to  this  program  as  well  as  that  in  Appendix  A.  The 
computational  program  consists  of  a  main  program  and  a  sub¬ 
routine.  A  photographically  reproduced  copy  of  the  listing 
for  these  programs  is  given  following  the  list  of  symbols. 


Program 
Symbols : 


TBS  Maximum  peak  temperature  in  degrees 

Centigrade . 

TP1  First  pulse  width  in  seconds. 

TEMPI  Maximum  peak  temperature  in  degrees 

Centigrade  which  corresponds  to  a 
pulse  width  of  TP1. 

TP2  Second  pulse  width  in  seconds. 

TEMP2  Maximum  peak  temperature  in  degrees 

Centigrade  which  corresponds  to  a 
pulse  width  of  TP2 . 


TP1  and  TEM1  may  both  be  zeros.  If  it  is  known  that  one  pair 
of  estimates  is  better  than  the  other  this  pair  should  be 
selected  as  TP2  and  TEM2. 
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PAIN  PROGRAM  FCR  FINDING  PULSE  LENGTH  WHICH  WOULD  GIVE 
FINAL  TEMPERATURE  USES  SUBROUTINE  TEMPITP.T2) 

COMMON  TMAX,PC,NMAX,GAM,ALF,CP,CUR,RZ,RATE 
5  RE AC  INPUT  TAPE  0,101,  TMAX , PC , NMA X , GAM , AL F , C P 
1C  R  E  AC  INPUT  TAPE  C  ,  1 C 2  ,  CUR , R Z , R A T E , TB S , T P2 , TEM2 , T P 1 , TE M 1 

11  IF(RATE)  40,40,12 

12  N  =  0 

2b  TP=TP2+(TBS-TEM2)«(TP2-TP1)/(TEM2-TEM1) 

30  T  P 1  =  TP2 

31  TEM 1=TEM2 

32  T  P  2=  TP 

33  N  =  N*1 

34  CALL  TEMP1TP2.TEM2) 

3b  IF(N-b)  38,36, 10 

38  IF (ABSF ( TEM2-TBS1-5.0)  10,10,25 

40  CALL  RETURN 
4b  STOP 

1C1  FORMAT (F8.C.F8 .4 , 18, 3E  1  1 .4  1 
102  FORMAT  (4F8.0,2(E 1 1.4.F8.2 ) ) 


C  SUBROUTINE  FCR  USE  WITH  PRECEDING  PROGRAM 

SUBROUTINE  TEMPITP.T2) 

COMMON  TMAX, PC»NNAX, GAM, ALF, CP, CUR.RZ, RATE 
7  NN  =  C 

12  IF (  1  •  0  -  T  P  »R  A  T  E  )  115,115,16 

16  T 1=0.0 

17  R  1  =  R  Z 

19  TLST  =  T  I 

20  N  =  0 

21  FNG=0. 0 

23  TS= 1.0/RATE- TP 

24  EXP=0.36787944«»(TS«GAM/CP1 

28  WRITE  OUTPUT  TAPE  0,304,  CUR , RZ , R A T E , TP , G AM , AL F , C P 

29  NN=NN*5 

3 C  C  SC  =  COR • »2 
31  ENZ=RZ»CSu 

40  EN=R1*CSC 

41  D=EN»ALF-GAM 

42  C=CP/D 

43  E  =  2. 71828 1E3«« ( TP/C  ) 

44  A=EN/D 

45  N- N ♦  1 

46  B  =  A • AL  F 
4e  E  C  =  E -  1 . 0 

50  T2  =  T 1 ♦ E  D • A 

51  T3=T2»EXP 

52  R 2  =  R Z • t  1 .0*ALF»T2  ) 

53  R3  =  R Z • (  1  .O+ALF  «T3  1 

60  PENG=EN«TP*B • l C «EC-TP  )  *ENZ 
62  ENG  =  ENG  +  PENG 

74  AN  =N 

75  T l ME  =  AN / R  A  T  E 

76  IF ( N  N—  4  0  )  79,79,77 

77  WRITE  OUTPUT  TAPE  0,307,  N.PENG 
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NN  =  C 

7 8  GC  TO  60 

79  WRITE  OUTPUT  TAPE  0,308,  N , PENG 

80  WRITE  OUTPUT  TAPE  0 , 309  , T I  PE , ENG 

81  WRITE  OUTPUT  TAPE  0,310 

82  WRITE  OUTPUT  TAPE  0,311,TJ,R1 

83  WRITE  OUTPUT  TAPE  0,312,T2,R2 

84  PX=(T2+273.C ) / ( TL S T+ 2 73 . 0 > 

86  IFITMAX-T2)  89,89,87 

87  IF(PC-PX)  86,89,89 

88  IFINMAX-N)  89,89,90 

89  RETURN 

90  R  1  =R 3 

91  TL  ST  =  T2 

92  T  1  =  T  3 

93  NN=NN+6 
95  GC  TO  40 

1  15  CALL  RETURN 
116  STOP 

304  FORMAT  l  1H164I  IH  )  1 8H  THE  RM  AL  P  A  R  AM  E  T  E  R  S/ 1  0  I  1 H  )  78HCURRENT  RESIST 
1ANCE  PULSES  PULSE  WIDTH  GAMMA  ALPHA  CEE  PEE/1211 

2H  ) 78HAMPS  OHMS  PER  SEC  MICROSECONDS  MICWTS/DEG  P/ 

3M/DEG  MICJLS/CEG/F 16.2.F 13. 3,F 1 1 .2,6PF 14.3.2F 1 1 .  1  ,F 1 2.3) 

307  FORMAT (  lH  I  16 (  1H  11IHCVCLE  COUN T  I  1 0 ,  1  3  I  1 H  ) 15HENERGV  IN  PULSE3PF11. 
13.12H  MILLIJCLLES) 

308  FORMAT  I  1H0  16 (  1H  111HCYCLE  COUNT  I  1  0  ,  1  3  I  1 H  )  15HENERGY  IN  PULSE3PF11. 
13, 12H  MILLIJCLLES) 

309  FORMA  T  I  1 7  I  1 H  )  12HELAPSEC  T I ME3 PF  1  2 . 2 , 25H  MILLISEC  TOTAL  ENERGY 
IF  I  1 .3, 12H  MILLIJCULES) 

310  FORMAT {  1 7 {  IH  )  1 1 H T EM PER A TUR E2 3  I  1 H  )  1  OHRE S I  STANCE ) 

311  FORMAT  I  191  IH  )  1 2HBEF0RE  PUL SE F 1 C . 2 , 24 H  DEG  C  BEFORE  PULSE 

1F12.3.5H  OHMS) 

312  FORMAT (  191  IH  1  1  lhAFTER  PUL SE F 1 1 . 2 , 2 3H  DEG  C  AFTER  PULSE 

1F13.3.5H  CHMS) 
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